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Abstract. A long memory process has self-similarity or scale- invariant 
properties in low frequencies. We prove that the log of the scale-dependent 
wavelet variance for a long memory process is asymptotically propor- 
tional to scales by using the Taylor expansion of wavelet variances. 
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1 Introduction 

Many natural phenomena have shown self-similarity or fractal properties in di- 
verse areas. It is well known that time series in neuroimaging such as electroen- 
cephalogram, magnetoencephalogram, and functional MRI have fractal prop- 
erties |1I2| . Most of fractal time series can be well modeled as long memory 
processes. 

A long memory process generally has scale-invariant properties in low fre- 
quencies; for example, the scale-dependent wavelet variance has asymptotically 
power law relationship in low-frequency scales, which implies its fractal property 
[1] . Such a property also has been demonstrated by Achard et al. [3] by exploit- 
ing the Taylor expansion of wavelet variances; regrettably, they did not provide 
the detailed proof. Here, we show the detailed procedure of Achard et al.'s proof. 

2 Wavelet Variance of Long Memory Process 

A long memory process can be denned as follows as suggested by Moulines et 
al. @]; that is, a real-valued discrete process X := {X(t)} t _ 1 ... N is said to have 
long memory if it is covariance stationary and have the spectral density 

Sx(/) = |l-e-^r 2d Si(/) (!) 

for < d < 1/2 where 5*x (/) is a non-negative symmetric function and is 
bounded on (— 7T, 7r). 
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Theorem 1. Let Wj t k is the wavelet coefficient of the time series X at the j-th 
scale and the k-th location, and let T4j{f) be the squared gain function of the 
wavelet filter such that 

2 j for 2tt/23 +1 < \ f\ < 2n/2^ 
erwise 



1 othe^"-"" 



Then, := Var(Wj^) satisfies 

j(j) = 2n % 3 {f)S^{f)df. (3) 

J — 7T 

Achard e£ a/. [5] made its Taylor expansion to estimate the long memory param- 
eter of X as shown in the corollary [TJ We give the detailed proof of the corollary 

CD 



Corollary 1. Let (/) = 1 + /3/ 2 + o (f) . The 

r 2ir/2 j 



7(i) = 2^ +1 J2l +1 Re W)) df = K2 2d > (l + a 2 l- + o (± 



2 J 6 V 2 

Let us define four integrals: 



(4) 



where 

1 - l/2 1 - 2d 9 ,J 
^ = 2 (27r) 2 - 2d , (5) 



r9 ,2 fd \ (1-1/2^/ (3 -2d) 
fl2 = (27r) V 12 +/3 J (1-1/21-") /(I -2d) ' (6) 

Proof. It is well known we have the following Taylor series 

sin- M (//2) = (//2)" 2d + f (f J +o(/ 3 ). (7) 
From (HJ and ©, 

^ (/) = 2" 2d ( (I) ~ M + 2 4 (?) 2il ~ d) ) (1 + /3/ 2 + o ( f)) . (8) 



2I/2 „,(0 ^= ■ ,10) 

2tt/2 j /y\2(l-d) ^3 (-22^ _ g~) 2 i( 2d -3) 

27T/2J + 1 



2; rf/ 4(2d-3)7r 2rf ' (U) 
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h ■= 



r 27r/2 J 



'2ir/2J-l 



f\ 2(l-d) 
J \ J-2 



it 5 (2 2d - 32) 2^ 2d - 5 ) / 
f = 4 (2d -5)^ :=^ (2 " 5) - (12) 



Since if h (Af) for Af ->■ then / a a+4/ & (/}/) df = o((4/) 2 ) 



2tt/2 



-2(/ 



_(/ a )d/ = 2 2 *o( 4- 



2tt/23 

^2^/2^ /,x2(l-d) 

/6:= / i 

From dHJ and the integrals ©-([HI), 

<.2tt/2 j ' 

1 (j)=2irV +1 Re(5(/))df 

J2n/2i+ 1 



5 (/ 2 )d/ = 2 2 *o 



2 4 J 



2d 



h + Ph + + 2 4$h + h + =£h 

6 6 6 



(13) 
(14) 



Note that I 3 = I 2 /4 and 



->2d 



2 2d - 2 2d - 3 2 2 i 
Thus, 

7 (j) = 7 r2^+ 2 ( 1 - d ) 
= 7 r2-'' +2 ( 1 - d ) 



l_ 2 2d-i 3 -2d 2 2 J 



(15) 



2 2 J'/ 



2-V 



(i + n) a^w-v + 2 A pA ^ 2d -v + 2 2 * r ^ 



where 



a 2 := 09 + d/12)Mi2 
= (2-) 2 (£ + /3 



a , d \ (1 - l/2 3 - 2d ) / (3 - 2d) 



12 



(l-l/2 1 - 2d )/(l-2d) 



Then, we can have 

7 (j) = i^2 2 * 
= K2 2 * 



1 2dA 4o 1 1/1 

1 + fl2 2W + y^w + Flw 
1 + fl2 w + (w 
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where 

1 — 1 /oi-2d 

K := 7rA2 2(1 " d) = 2 7 , (27r) 2 - 2rf . (16) 
1 - 2d y ' v 1 

In the corollary ([T]), as the scale j increases, the wavelet variance con- 
verges; in other words, 

j(j)^K2 2 <V. (17) 

It implies that we can simply estimate the long memory parameter d by linear 
regression method if we have the estimates of wavelet variances. 

3 Discussion 

We verified Achard et al.'a proof on asymptotic properties of wavelet vari- 
ances of a long memory process. Unfortunately, Achard et al.'s Taylor expan- 
sion of wavelet variance is not perfectly consistent with our result; while K = 

2 1 ~\^ 2 d — (27r) 2 ~ 2d in our result, they computed as K = 2 1 ~ 1 { 2 2rf — (27r) 1_2d . 

On the other hand, the assumption on short memory such that (/) = 
1 + Pf 2 + (/ 2 ) loses generality of their proof since there exist more general 
classes of short memory; indeed, Moulines et al. proved that the asymptotics of 
wavelet variances hold when the short memory (f) belongs to the function 
set 7i {fi,L) defined as the set of even non-negative functions g on [—1/2, 1/2] 
such that 

\g(f)-g(0)\<Lg(0)\27rff. (18) 

Nevertheless, their method is relatively simpler than other mathematical proofs. 
Moreover, we will be able to apply their method to investigate the asymptotic 
properties of multivariate long memory processes as Achard et al. already at- 
tempted for bivariate long memory processes [3|. 

References 

1. E. Bullmore, J. Fadili, V. Maxim, L. Sendur, B. Whitcher, J. Suckling, et al., 
Wavelets and functional magnetic resonance imaging of the human brain., Neu- 
rolmage. 23 Suppl 1 (2004) S234-49. 

2. V. Maxim, L. §endur, J. Fadili, J. Suckling, R. Gould, R. Howard, et al., Fractional 
Gaussian noise, functional MRI and Alzheimer's disease, Neurolmage. 25 (2005) 
141-158. 

3. S. Achard, D.S. Bassett, A. Meyer-Lindenberg, E. Bullmore, Fractal connectivity 
of long-memory networks, Physical Review E. 77 (2008) 1-12. 

4. E. Moulines, F. Roueff, M. Taqqu, On the Spectral Density of the Wavelet Coef- 
ficients of Long-Memory Time Series with Application to the Log-Regression Es- 
timation of the Memory Parameter, Journal Of Time Series Analysis. 28 (2007) 
155-187. 

5. D.B. Percival, AT. Walden, Wavelet methods for time series analysis, Cambridge 
University Press, 2000. 



